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Abstract 



We present here a calculation of the inertial mass of a moving vortex in 
cuprate superconductors. This is a poorly known basic quantity of obvious 
interest in vortex dynamics. The motion of a vortex causes a dipolar den- 
sity distortion and an associated electric field which is screened. The energy 
cost of the density distortion as well as the related screened electric field con- 
tribute to the vortex mass, which is small because of efficient screening. As 
a preliminary, we present a discussion and calculation of the vortex mass us- 
ing a microscopically derivable phase-only action functional for the far region 
which shows that the contribution from the far region is negligible, and that 
most of it arises from the (small) core region of the vortex. A calculation 
based on a phenomenological Ginzburg-Landau functional is performed in 
the core region. Unfortunately such a calculation is unreliable, the reasons 
for it are discussed. A credible calculation of the vortex mass thus requires 
a fully microscopic, non-coarse grained theory. This is developed, and re- 
sults are presented for a s-wave BCS like gap, with parameters appropriate 
to the cuprates. The mass, about 0.5 m e per layer, for magnetic field along 
the c axis, arises from deformation of quasiparticle states bound in the core, 
and screening effects mentioned above. We discuss earlier results, possible 
extensions to d-wave symmetry, and observability of effects dependent on the 
inertial mass. 
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1. Introduction 

The discovery of high temperature superconductors has led to a renewed interest in the 
mixed phase. Several novel phenomena arising from their short coherence length, layered 
nature and large superconducting transition temperatures have been theoretically and ex- 
perimentally studied(l). One area of interest is the existence of effects connected with vortex 
dynamics, e.g. quantum creep (2,3), anomalies in the Hall effect (4-6), and ac electromag- 
netic response (7,8). These phenomena are not fully understood, partly because of the lack 
of a well developed first principles theory of vortex dynamics, especially in the quantum, 
interacting vortex regime. A number of recent contributions address parts of the problem 
(9-12), especially the Magnus force driven dynamics in the presence of dissipation. 

A necessary ingredient in all considerations of the motion of a vortex is its inertial mass. 
This quantity, generally believed to be small, is surprisingly ill known and its origin is not 
well understood. (See Ref. 1, for example). Not much attention has been paid to this 
question because, for some phenomena, the dynamics is governed by the large dissipation 
(13) or the strong Magnus force (9-12) and the inertial mass could be irrelevant. However, 
there is experimental evidence for a low dissipation regime in cuprate superconductors (5), 
and for a Magnus force smaller (14) than standard estimates (9-12). It is thus quite possible 
that the inertial mass could affect dynamical processes involving vortices. (These questions 
are taken up in Section IV). Also, in the absence of an understanding of what contributes 
to the vortex mass, and how much, it is difficult to meaningfully discuss the question of 
whether or how such a mass influences vortex dynamics. We therefore present here an 
extensive discussion and a calculation of the vortex mass, and go into the question of mass 
related phenomena in Section IV. This work was first reported in 1994(15). 

We first estimate the mass using a phase-only functional (Section II) which is known to 
give a good description of the system far from the vortex core, where the amplitude of the 
order parameter is nearly constant and the only relevant degree of freedom is the phase. 
Recently, Duan and Leggett (17) and Duan (18) have given a careful discussion of this 
approach where the time dependent order parameter (phase) of a moving vortex causes the 
electronic density to fluctuate. This in turn gives rise to an electric field which is screened. 
The energy of the density distortion and electric field energy are the cause of the mass. The 
phase-only functional used has been derived microscopically (19) and the result obtained 
thus has a microscopic significance and gives an accurate estimate of the contribution to the 
vortex mass which accrues from transitions induced in the electronic scattering states of a 
vortex by the vortex motion. The contribution of this process from the far region turns out 
to be very small due to the efficient screening; thus most of the mass comes from the core 
of the vortex. 

We then calculate the core contribution to the mass using a phenomenological Ginzburg- 
Landau functional as has been conventional since the early work of Suhl (16). This functional 
is not derivable microscopically and is used mainly as an interpolating formula which reduces 
to the correct phase-only functional in the far region. However, obviously the coarse grained 
GL approach is unrealistic for effects within the core which is of the same size (£) as the 
coarse graining scale of the theory. Further the screening length in GL theory is proportional 
to where ip is the superconducting order parameter, and thus diverges at the core. 

This is clearly an artefact of the GL approach, as the efficiency of Coulomb screening is 
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related to the electronic compressibility, and is expected to be largely independent of the 
superconducting order parameter. The GL estimate of the vortex mass is clearly unreliable 
and is presented here mainly to contrast the correct microscopic calculation of the mass 
which forms the main body of this paper. The correct, microscopically obtained vortex 
mass, is very different from the GL estimate in its dependence on the basic parameters 
of a superconductor. Inspite of this, the numerical value of the mass obtained in the GL 
approach is (largely by accident) in the same range as that obtained from the microscopic 
theory for parameters appropriate to the cuprates. 

We then present the first correct microscopic calculation of the vortex mass (Section 
III). We employ the self consistent pair field approximation which has been used extensively 
for static vortex structure, quasiparticle energy levels, etc. (20-23). We make a Galilean 
transformation to a frame of reference where the vortex is at rest. In this frame, the motion 
of the vortex acts like a perturbation of the form u.fP p where u is the vortex velocity and fP p 
is the momentum operator for the electrons. The inertial mass is obtained by integrating out 
the electronic degrees of freedom to second order in u. The coefficient of the (u 2 /2) term in 
the effective action is the effective mass of the vortex. The mass is found to originate from a 
polarization process involving the virtual excitation of the lowest energy quasiparticles in the 
bound electronic states (quasiparticle states) localized in the core of the vortex. The small 
core size (£ ~ 15 A) in the cuprate superconductors implies that the lowest unoccupied state 
is separated by a sizeable gap (~ 100K) from the highest occupied state below the Fermi 
level (23), in strong contrast to conventional superconductors. The existence of this large gap 
in the core quasiparticle spectrum has been recently observed (43) by STM measurements 
in the vortex core region in Y Ba2Cu2,0-js- This virtual transition process gives rise to a 
large vortex mass (m* ~ 25m e ). However, strong dielectric screening drastically reduces the 
mass and leads to a value m* ~ 0.5m e per Cu0 2 layer. We discuss the physical reason for 
a mass of this size in terms of the basic length scales and the screening process. 

In the final section (Section IV) we discuss the calculation critically, compare with other 
results, consider the calculation of a vortex mass for a non s-wave superconductor and go 
into the question of when effects due to the small vortex mass might be observable. 

II. Ginzburg Landau calculation of the vortex effective mass 

The most natural way of discussing the motion of a singularity in the phase 8 of the 
superconducting order parameter is the time dependent Ginzburg Landau theory (17 - 19) 
where the free energy (or action) is expressed as a functional of the phase of the supercon- 
ducting order parameter. This functional provides a good description of the region far from 
the centre of the vortex where the amplitude of the superconducting gap is nearly constant. 
We describe the functional and briefly summarize known results for the mass contribution 
from the region outside the core (17,18). 

Outside the core, a phase only Hamiltonian is sufficient, and the action functional S per 
length L (at T = 0) is given by 

{S/L) = S e + S em (la) 



where 

/+oo r o„ 4 \ 2 

dt I df 
-oo 



h ) 2 \ he h 



(16) 
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and 

r (Viof-(VxI) ; 



/+oo r 
dt / dr 
-oo J 



8tt 



(lc) 



where r is the coordinate in the plane perpendicular to the magnetic field. In Eq. (lb), the 
first two terms are the energies associated with pair charge and pair current (or velocity) 
fluctuations respectively. The term linear in 9 is a total time derivative which has no physical 
consequences in the absence of vortices. Its physical origin can be understood (10) in terms 
of the Berry phase associated with the adiabatic motion of a vortex and it gives rise to the 
Magnus force. Notice the absence of any term linear in Aq which is constrained to be zero 
because of charge neutrality in the electron-ion system. 

The action functional of Eq. (la) can be obtained microscopically by starting with 
electrons interacting with a pair potential (whose magnitude is nearly constant in the far 
region), going to a gauge where the order parameter is real and then integrating out the 
fermions. The details of the derivation have been outlined in (19). In Fig. (1), we show the 
Feynman diagrams which contribute to the action of Eq. (la) at T = in the clean limit. 
The coefficients a,\ and a 2 are the appropriate polarizabilities. For a weakly interacting, 
clean Fermi gas, they have the values 

(2a) 




and 
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where X^ F (the Thomas Fermi screening length) is given by \*rfp = (Qime 2 /e^) and \ L (the 
London penetration depth) is given by A^ 2 = {Anne 2 ) /{m e c 2 ). A and A are the scalar 
and vector electromagnetic potentials respectively, the two terms in Eq. (lc) being just the 
field energies. In case there is an external magnetic field Hq, the term |(V x A) 2 /87r| is 
modified to {(V x A — H ) 2 /8ir}. The functional is clearly gauge invariant. In Eq. (1) the 
displacement current term {\^) 2 has been omitted because of the largeness of c. 

The motion of the vortex in a charged superconductor gives rise to an electric field (or 
potential A ). To find this A , we minimize S with respect to A which yields 

v% -2e ai u_2jM (3) 



47r h V h 
Now for a vortex moving with a (small) uniform velocity u, we make the assumption that 

9(r,t) =9 (r-ut) (4) 

where 9 (f) is the phase around a static vortex at the origin. Thus we have 

0(f, t) = -u.V0o (5) 

Using this in combination with the known 9q, one finds the potential Aq from Eq. (3) and 
the extra energy from Eq. (1). The result is 

,_, Train 2 , , £ _2 — R~ 2 . ,„ . 

Eke = —hWl + ^-2 + fl _ 2 ) M 
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where R c is a long distance cutoff (whose magnitude will be taken to infinity in the end). 
Then it is easy to see that in the limit e — > which implies \ TF — > oo, the expression for 
Eke reduces to 7Ta ^ u In(^) which diverges logarithmically as the long-distance cutoff R c is 
taken to oo. This is the expected result for a neutral superfluid (17, 18). For the present 
case of a charged superconductor we are always in the limit \ TF <C £ and the expression for 
the kinetic energy simplifies to 



J KE 



Traiii 2 \ FF 

The expression in Eq. (6b) can be rewritten as 




(66) 




Eke = E sv \^\\^\\—^-\ (6c) 



where E sv = (^^) 2 ln(^j-) is the London energy (per unit length) of a static vortex which 
results from the terms involving transverse fluctuations in the current and magnetic field in 
the phase functional of Eq. (1). The result is then easily understood on physical grounds; 
the velocity u is to be compared with the natural velocity scale v F of the Fermi system. 
The second bracketed factor in Eq. (6c) is due to the reduction of charge fluctuations by 
screening. The former occur on a length scale £, while the screening length is X TF <C £. On 
substituting appropriate numbers in Eq. (6b), i.e. £ ~ 15A and \ FF ~ 1A we find that 
the contribution of this source to m*j ~ 7 x 10~ 4 m e per layer, (m e being the electron mass) 
assuming an interlayer spacing d ~ 10A This is an extremely small number; one reason for 
its smallness is the screening factor (A^ F /£) 2 ~ (1/200). The result, (Eq. (6b)), has been 
obtained earlier by Duan (18), where details may be found. 

The functional Sg of Eq. (1) is not Galilean invariant and there are additional terms 
involving higher order derivatives of 9, whose inclusion is necessary to restore this sym- 
metry. These arise from the fact that the Gaililean invariant combination involving the 
time dependent order parameter phase 9 is the local electrochemical potential 8[j,(r,t) = 

/ HV8 eA\2 

— 7}9 + eA — 2 2m c The correct functional was recently obtained by Aitchison and co- 
workers (24) and is of the form 



/+oo r 
dt / df 
-oo J 



*[(6 - ^) + MW ~tf )2 f - ^9 + (V9 - 2 -^f] 



(7a) 



The additional diagrams contributing here, are shown in Fig. (2). The harmonic electric 
field terms in Eqs. (7a) and (lc) can be integrated out to give the action functional 



It is clear from Eq. (7b), that additional contributions accruing to the vortex mass come from 
the coupling of charge fluctuations to the supercurrent fluctuations leading to non-adiabatic 
corrections in the supercurrent distribution that are proportional to the vortex velocity. 
However, an explicit calculation shows that all such contributions to the vortex mass are 
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screened at least doubly more efficiently and are therefore smaller than that estimated above 
(Eq. (6b)) by a factor of (Atf/£) 2 ~ -005 . Since the far mass estimated before is already 
small, this additional contribution is smaller still, and therefore negligible. In fact in the 
perfect screening approximation, where the total electro-chemical potential 5fi is set equal 
to zero locally, we get back the result of Eq. (6), the contribution of the additional terms 
being zero as expected. It is worth emphasizing that the results obtained from the phase- 
only functional are essentially microscopic (see the Feynman diagrams in Fig. 1 and Fig. 
2) and give an accurate estimate of the contribution to the vortex mass from the electronic 
scattering states which are extended in nature and live mainly outside the core. 

We now consider the contribution to the vortex mass arising from the core. The ampli- 
tude of the superconducting order parameter changes as a function of radial distance from 
the center of the vortex in the core region. We use a phenomenological Ginzburg - Landau 
action functional per unit length 



dt \ dr 



2mv F 



'hd_ 



-2eA )1>(r,t) 



1 

2m 



2e T \ , 

-V A i\) 

% c J 



+ 



dr 



(VA ) 2 - (V x Af 
8tt 



(8) 



to estimate the core contribution to the vortex mass. 

This functional is not derivable microscopically and has been chosen primarily as an 
interpolating formula which reduces to the correct phase-only functional in the far region in 
the limit ^ — > ^Jr^ s e %e and gives the correct (linear in r) dependence for the amplitude of 
the order parameter near the center of the vortex. The parameter m which appears above 
is therefore fixed by requiring this functional to reduce to the action functional of Eq. (1) 
in the "phase-only" approximation. For the potential V(|-?/>| 2 ) we assume the standard form 

vW) = !H 2 + fH 4 - 

Starting with the early work of Suhl (16) and subsequent work by others (see Ref. (1) 
and references therein) all estimates of the vortex mass have proceeded from this functional. 
The value of the mass obtained from this functional is not expected to be very accurate 
for reasons that are discussed below. However we present a calculation of the vortex mass 
from this functional mainly to contrast and highlight the microscopic calculation presented 
in Sec. III. 

We again (see Eq. (4)) make the ansatz that for small velocity u, the vortex moves 
rigidly i.e. 

*P(r,t)=Mr-ut) (9) 

where ipoif) is the order parameter configuration associated with a static vortex. We further 
assume a common and fairly accurate explicit form for ipo(r), namely 



^aif) — v^-s tanh(r/£) exp(i(, 



(10) 

where n s is the superfluid density far from the vortex, and (r, <fi) are the radial and angular 
coordinates of the two dimensional vector r. Again, minimizing Eq. (7) with respect to A , 
we find 

V 2 A 3eh ( , dil>* ,*dV>\ 3 



47T 



imvp 



dt 




mvp 



(2e) 2 \^\ 2 A 



(11) 
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. Using this condition in the action functional of Eq. (8) we find that the extra energy due 
to vortex motion is 



AE 



I 



dr 



3ft 2 
2mv 2 F 



dip 
W 



3ehA 
2mv\i 



(12). 



The first term in Eq. (12) is due to the density fluctuations induced by vortex motion, 
and the second describes the reduction due to screening. The first term is readily computed 
using the ansatz Eq. (9) for the time dependence of ip(r,t) and the form Eq. (10) for the 
coordinate dependence of ipo(r). We are interested here only in the core contribution to the 
mass, the contribution from the far region having been previously determined (Eq. (6)). 
The radial integration in Eq. (12) is therefore performed over the range < r < £ (core 
region). It gives a mass per unit length 



m unscreened 



- m u • = 0.61 ( — 



(13) 



where a = is the Bohr radius, and (X TF ) 2 = 4irn s (2e) 2 / (mv F /3). For the cuprate 



superconductors, we find that the unscreened core vortex mass per layer is 



TO, 



0.19m P 



(14) 



An unexpected, and incorrect, feature of Eq. (13) is the lack of dependence of the vortex 
mass on the core size £. This is a consequence of the fact that the core mass is proportional 
to the gradient energy per unit length in the core, which is scale invariant and does not 
depend on the natural length scale £ in two dimensions. This is inevitable in any local 
continuum free energy theory. 

We now consider the second or screening term in Eq. (12). One clearly needs to know 
the electric potential A induced by vortex motion. Setting 



Ao(r) = v{r)u • 



(15) 



we find that v(r) satisfies the radial equation 



d 2 



^19 1 



\if>(r)\' 



dr 2 r dr 



(A 



TF 



) 2 



v(r) 



l_h_ 

4vr2e 



A 



TF 



tanh 2 (r/£) 



(16) 



This is just the radial Poisson equation with a screening term (last term on the left hand 
side), and a source term (on the right hand side) appropriate for a two dimensional system. 
The 'effective screening length' XtfV^s/I^^I ^ s r dependent, and increases as \ip(r)\ de- 
creases, i.e. with r — > 0. This is yet another unrealistic feature of the Ginzburg Landau 
theory, since one expects screening which is related to the electronic compressibility, to be 
relatively independent of superconducting order. Eq. (16) is solved numerically, with ap- 
propriate boundary conditions at r = 00 and r = 0. The v(r) and the A (r) thus obtained 
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(see Eq. (15)) are used to calculate the second term in the vortex core energy (Eq. (12)). 
Adding the contributions of both the terms, the final result for the mass per layer is 



m* (core) = m] 



Hii;H3 1061 - 041 - 



(17) 



where d is the interlayer separation. In Eq. (17), the first term in the square brackets is due 
to the density distortion in the core, and the second term is the negligibly small correction 
due to screening. Thus the total mass per layer, for a cuprate superconductor at T = OK is 



In the absence of a microscopic theory, the phenomenological Ginzburg Landau approach 
has been used, mainly as a dimensional aid, to estimate the vortex inertial mass. The values 
quoted ( Ref. (1) ) are in the range 0.2 to 2 m e , and turn out to be (largely by accident) 
not far from our microscopic result (see below) for parameters appropriate to the cuprates. 
As a preliminary to the microscopic calculation which forms the main result of this paper, 
we have performed a detailed calculation using this functional to bring out its inadequacies. 
The above detailed analysis of the different contributions to the mass shows that within 
the phenomenological GL theory it is due to the essentially unscreened density distortion 
induced in the core. The unscreened mass estimated above too is incorrect as the gradient 
expansion implicit in a phenomenological theory like the GL functional of Eq. (8), breaks 
down at short length scales and strongly non-local (in space) effects lead to a much larger 
mass (see Sec. III. below) in the unscreened case than calculated here. This picture is 
clearly incorrect on another count as well. Electronic screening processes are not expected 
to be affected much by the onset of superfluid order and a strong reduction of the mass is 
expected because of Coulomb screening. This aspect too is explicitly seen in the microscopic 
theory. In short, the microscopic calculation shows that the phenomenological GL picture 
is wrong on all counts, as we shall see now. 
III. Microscopic calculation of the vortex mass 

We present now a microscopic calculation of the vortex inertial mass for a layered su- 
perconductor at T = 0. The dynamics of the system of paired electrons is described by the 
action 



m* ~ 0.19m e 



(18). 



S 




em 



(19a) 



where 



L{ (f, t) = -eA (f, Id) J2 ^a(f, WU?, t) 




(196) 



(19c) 



a 



L? air (r,t) = -lMr,t)Mr,t)Mr,t) + h.c] - 



Mm 2 



(19d) 



V 
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and 




df 



jyAp) 2 - (V x A) 2 
r 8tt 



(19e) 



The electrons at (f, t) on layer /, with spin a are represented by the Grassmann field 
variables 4>i a (f, t), ipi a (f, t). The electronic kinetic energy is described by the term h[ (r, t), 
i.e. Eq. (19b) and the electronic coupling to the Coulomb potential by the term L{(f,t) 
(Eq. (19c)). Eq. (19d) is the mean pair field decomposition of the Cooper pair attraction, 
appropriate for a BCS s-wave superconductor and V is the strength of the attractive contact 
interaction. The last term S em (Eq. (19e)) is the electro-magnetic field contribution. 

Several possible modifications and generalizations, such as order parameter symmetries 
other than s-wave and effects due to quasiparticles at T 7^ are briefly discussed in Section 
IV. We have neglected above the small interlayer (Josephson) coupling, so that the (pancake) 
vortices in different layers 

are coupled only via electric and magnetic fields. This neglect has very little effect on 
the vortex mass, which is overwhelmingly due to processes occuring within a layer, and as 
a matter of fact, to processes within the small core. 

Every layer has one pancake vortex. All the pancake vortices are assumed to lie along a 
straight line parallel to the magnetic field which is perpendicular to the layers and to move 
with a uniform velocity u, so that the vortex (core) coordinate is (Ro,ld) where R = ut is 
a vector in a plane parallel to the layers and the z co-ordinate Id specifies the position of 
the layer. The pair potential A;(r, t) in the presence of such a uniformly moving vortex is 
a function of (f — R (t)). It can thus be written in an adiabatic approximation as 



Here Ao(r) is the magnitude of the pair potential and 0(f) the phase, for a static vortex 
situated at f — 0. The pair potential does not depend on I. The action functional relevant for 
vortex dynamics is found by an expansion of the microscopic action in powers of the vortex 
velocity u after integrating out the electronic degrees of freedom. The vortex mass is then 
determined from the term quadratic in u (the dissipative and Magnus forces coming from the 
linear term). This calculation is most conveniently carried out in the rest frame of the vortex. 
In this frame, the pair potential seen by the electrons is Ao(r) exp[i(0(f) — 2m ^' r + ^-^)] 
i.e. the pair potential corresponding to a vortex at rest with extra phase factors coming 
from the macroscopic supercurrent and kinetic energy of the electrons which acquire an 
extra velocity — u in this frame. The Coulomb potential seen by the electrons in this frame 
becomes A — ^ • A while the vector potential is the same as before. All physical quantities 
can be calculated in this frame and then transformed to the lab frame as necessary. Note 
that this does not assume Galilean invariance. 

After a gauge transformation, the action for the system in this moving frame can be 
written as 



A t (f,t) = A (|f - R (t)\) e ^-^W) 



(20). 



S = So + S 1 



(21) 



where the static vortex action is 



S = dt dr^ L{(f,t) - A (r)[e ie ^(r,t)iP n (r,t) + h.c] 




V 
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— J dt J dz J dr- 




(22a) 



and the perturbation due to vortex motion is contained in 




dr 



,(VA,) 2 



(226) 



In Eq. (22b), pi(f,t) is the momentum density operator for the /th layer, 



Pi(r,t) 



J2i>ia(?,t)-Vi> la (r,t) 



(23a) 



a 



and pi(r,t) is the density operator 



Pi(r,t) 



^24> la (r,t)ij{ a (r,t) 



(236) 



The first term in Eq. (22b) is linear in u, coupling to the electron momentum. The second 
term is the electric potential energy of the (nonuniform) electron density around the moving 
vortex. The associated electric potential has to be determined self consistently. 

Since the London screening length (A) is much larger than the core size (£), the magnetic 
field in the core is nearly the same as the external magnetic field, deviations from this being 
of order (£/A) 2 . However, the vector potential associated with this field is negligible in 
comparison with the gradient of the phase of the superconducting order parameter and has 
been ignored in the following calculations. The former « Hr whereas the latter <j> /r. 
Thus for the core region (r < £), we find that for H <C 0o/£ 2 the vector potential may be 
ignored (20). Corrections due to the vector potential can be estimated to be of order H/ H c2 
where H is the external magnetic field and H c2 is the upper critical field and are therefore 
small in the dilute vortex limit. 

Now a systematic expansion in powers of u becomes possible by integrating out the 
electrons giving rise to the Lagrangian that describes vortex dynamics. Dynamics of vortices 
(classical or quantum) can be studied either by working directly with the vortex action or 
alternately by introducing a canonically conjugate momentum which permits one to go over 
to the Hamiltonian formalism. 

The Magnus and dissipative forces come from the term linear in u and can be obtained 
by taking the gradient of this term with respect to the vortex co-ordinate (37). The inertial 
term in the action is quadratic in the vortex velocity. We are thus interested in the change of 
action to second order in u, the coefficient of {u 2 /2) in the change being the vortex effective 
mass to*. Clearly, this is calculable by going to second order in u and A . For ease of 
presentation and also to emphasize the importance of Coulomb screening in the core, we do 
this in two stages. First, we find the unscreened vortex mass, i.e. in Si (Eq. (22 b)), we turn 
off the Coulomb interaction by putting e = 0, and calculate the second order shift. We then 
calculate the effect of Coulomb interactions, i.e. the effect of the electric potential due to the 
electron charge density change consequent on vortex motion. We would like to emphasize 
that the unscreened mass calculated here corresponds to the contribution of the bound states 
localised in the vortex core. The contribution of bound states is finite. For a (hypothetical) 
neutral superconductor, this term is to be added to the contribution from the region outside 
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the core. This can be calculated from the (microscopically derived) phase-only functional 
(Eq. 1) by setting e = 0, and as mentioned earlier, is logarithmically divergent because 
longitudinal density fluctuations associated with vortex motion are not screened. Thus for 
the uncharged superconductor, the finite core contribution can be neglected in comparision 
with the log diverging far contribution. 

For a charged superconductor, the core contribution has to be calculated with the inclu- 
sion of screening effects which, we show ( Eq. (43) below) reduces the 'unscreened' or e = 
mass by a factor of fifty or more. To this we have to add the m* due to the far region, which 
because of Coulomb screening is finite, and is actually negligibly small (Eq. (6c)). 

The second order effective Lagrangian due to just the u.p term in Si is given using 
standard many body perturbation theory by ^-u 2 where 

m* = i J dr J dr' J dt (T [pf (r, t)pf(P, 0)] ) = - J dr J dr'xTi^ r ~*') ( 24 ) 

In Eq. (24), the vortex velocity has been assumed to be in the x-direction. Note that intra- 
layer averages such as Eq. (24) do not depend on the layer index. The correlation function in 
Eq. (24) is calculated with respect to the unperturbed action So (Eq. (22a)) which describes 
a single static vortex. This last problem of a static vortex has been studied extensively using 
the Bogoliubov-de Gennes self-consistent field theory (20-23). The electronic eigenstates, in 
the presence of a static vortex, are bound states which are localised and have appreciable 
amplitude only in the core of the vortex, or extended states which are scattered by the su- 
perfluid velocity and are primarily in the region outside the vortex core where the amplitude 
of the order parameter is nearly constant. The latter are well described by neglecting the 
variation in the amplitude of the order parameter and the scattering processes contributing 
to the mass involving these states are just those considered in the Feynman diagrams (Fig. 
1) contributing to the phase-only functional in Eq. (1). Thus the mass contribution from the 
deformation of the scattering states is accurately estimated by the calculation proceeding 
from the phase-only functional of Eq. (1) and has been shown to be negligibly small due to 
efficient screening. We will therefore concentrate in the following only on the contribution 
of the localised states to the correlation function of Eq. (24) to find the core contribution 
to the vortex mass. 

The eigenfunctions of the Bogoliubov-de Gennes equations for a single layer which are 
localised in the vortex core are the amplitudes u^r) and Vp(r), labelled by the (azimuthal) 
angular momentum quantum number fi because of the cylindrical symmetry of the single 
vortex problem. The Bogulibov amplitudes are related to the fermion field operators by the 
relations 

lh(r,f)= EMWt) (25a) 
alU 

and 

^\(r,t) = E^(07m(0 (256) 
all^t 

Here \i runs over all half-odd integers (positive as well as negative). The quasi-particle 
annihilation operators 7 M correspond to the empty (particle) states for ji > and filled 
(hole) states for /i < 0. In terms of the quasi-particle operators 7 M , the Hamiltonian for the 
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static vortex system described by the action So can be written in a diagonal form as, 



H = J2 e ^l/j. + constant (26a) 
M 

Our definition of the quasi-particle operators (Eq. (25)) differs from the standard definition 
(20) by a particle-hole transformation. In doing this, use has been made of the fact that 
if (u^r), Vfj,(r)) is an eigenstate with eigenvalue e M , the state (— v *(f), w* (f)) is also an 
eigenstate with eigenvalue — e M . Thus, 

(u,(f),v,(f)) = (-v*_^f),u*_^r)) (266) 

with 

e M = -e- M (26c) 

for « < 0. The ground state of the system has all states with « < occupied and all states 
with /i > empty. 

The correlation function of Eq. (24) is easily evaluated in this representation. The vortex 
mass can then be understood as arising from a polarization process involving a virtual (quasi- 
particle) transition from the highest occupied to the lowest unoccupied state. This process 
can also be viewed as a deformation of the ground state by the perturbation which mixes in 
higher energy states. In terms of u and v, the momentum-momentum correlation function 
can be written as 

i J dt(T[ P n?,t) P ?(p,o) 



II I) occ. occ. 



2 E E v) 



-1 



X 



M>0 /x'<0 



(27). 



The various terms here correspond to the different quasi-particle processes that contribute to 
the polarization. There is a simple selection rule for nonvanishing matrix elements, namely 
H = —fj,' = 1/2. We discuss this now. The details are worked out in Appendix I. The 
operator p x has an angular dependence of the form cos (f> (radial derivative term) and sin cf> 
(angular derivative term) where is the angle in the 2d plane of the vector r with respect 
to the x axis. The amplitudes u^(r), v^{f) can be written (20-23) as 

where f^(r) are functions only of the radial coordinate r. The dependence of u^r), v^(f) 
above implies that on integration of the matrix elements in Eq. (27) over the angle 0, the 
only nonvanishing terms are (//—//') = ±1. We also need one of the states /i to be unoccupied 
and the other // to be occupied. The only possibility among bound states is /i — — «' = \. To 
understand this, we exhibit in Fig. 3, the spectrum of eigenstates for parameters appropriate 
to the cuprate superconductor (23). The bound states within the gap have quantum numbers 
\i = ±|, ±|, The occupied (unoccupied) states have negative (positive) /i. It is now 
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clear that the only virtual transitions that satisfy the selection rule A/x = ±1, and are 
between occupied and unoccupied states, correspond to fj, — — // = |, i.e. the highest 
occupied and lowest unoccupied bound states. With this simplification used in Eq. (27) and 
the latter substituted into the expression Eq. (24) for tuq, we have 



41a I 2 

m* = (29a) 



ei/2 



where 

[ hdu 1/2 (r) ( s 

9x = J drv 1/2 (r)- — — (296) 

We thus need to know the energy 6\/ 2 and the core bound state wavefunctions fy 2 (r) in 
order to find m^. These have been determined self-consistently and numerically by Zhu et. 
al. (22). We use here the variational forms 

A„(r) = A tanh(r/0 (30a) 

f[ /2 (r) = A 1/2 J (k F r)e- r ^ (306) 

and 

f+ 2 {r)=A 1/2 J 1 {k F r)e- r ^ (30c) 

where 

V 2 = /*K(M0W?(M')]^ ! <3Q0 

is the normalisation factor. With £ = 15A, kp 1 = 3.36A and A = 60 meV, the expectation 
value of the energy < ei/ 2 > is 69 K, close to the self-consistent numerical value of 66 K 
obtained by Zhu et. al. (22). The wavefunctions are also very close. Using these, the 
expression ml can be evaluated (see Appendix I for details), giving a value 

ml ~ 25m e (31) 

The vortex mass obtained above can be estimated by the following simple physical argu- 
ments. The correlation function in Eq. (24) can be estimated as follows. Each of the 
momentum operators gives a factor of hkp which is the typical electronic momentum. The 
energy denominator of the correlation function (see Eq. (27)) is twice the bound state energy 
ei/2 which is the energy cost of the polarisation process involving creation of a 'particle-hole' 
pair. Finally there is a factor of two corresponding to electron spin. The bound state energy 

A 2 

(20) ei/2 is about The core contribution to the unscreened vortex mass (for parameters 
appropriate to the cuprates) is thus estimated to be 

h 2 k 2 

m* ~ £ ~ 100m e (32). 

This is larger than the value calculated above for the vortex mass (Eq. 31) by a factor 
of four. This discrepancy arises because the matrix element in the detailed calculation is 
smaller (by a factor of half approximately) than the dimensional estimate. 

We now consider the effect of Coulomb interactions. The dipolar charge distribution 
induced by vortex motion is screened efficiently by the electrons. This greatly reduces 
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the vortex kinetic energy. The reduction of the vortex kinetic energy is calculated using 
standard self-consistent linear response theory. The dipolar charge distribution produces 
an extra electric potential i.e. changes the electrochemical potential of the system. Any 
change in the electrochemical potential causes a change in density which therefore needs to 
be calculated self-consistently. In the following, we outline the calculation of the reduction 
of the vortex mass because of Coulomb screening. The details are provided in Appendix II. 
Varying the action in Eqs. (21) and (22) with respect to the electric potential Aq, we find 
at the extremum, the Poisson equation 



V 2 A (r, z) = -Ane ]T ( Pl (f)} 5(z - Id) 



(33). 



We find < pi(f) >, the electron density, to linear order in u and A to be 
(pi(r)) = e J dr'x°°(f, r')A (r',ld) - u j dr'-)fi x (r,r') 



(34). 



In Eq. (34), the density-density and density-current response functions x°° and X° x are 
given by 

(35) 



X °°(f,r') = -iJ dte(t)([ Pl (f,t), Pl (r',0) 



and 



i dte(t)( Pl (f,t), P f(?,o) 



(36) 



The second term on the right hand side of Eq. (34) is the source term for < pi >; it has 
to be determined by integrating the correlation function x^ x {^J r ') over the co-ordinate r'. 
This is evaluated to be 

' xT{ry)dr' = r 1 {r)\ (37) 



where 



and 



rj{r) 



A = 



\ e V2 



/l/2( r )/l/ 2 (0 Sm 



\ e V2 



h / dfv* 1/2 (r}- 



dx 



(38a) 



(386) 



It is thus clear from Eq. (37) substituted into Eq. (34) that the charge distribution generated 
by the vortex motion is proportional to sin0 i.e. it is dipolar in nature. We thus seek a 
self-consistent solution to the Poisson equation of the form 



A (r, z) = V(r, z) sin 



(39) 



Substituting this form in Eqs. (33) and (34), and using Eq. (37) we solve for A (r, z) (see 
Appendix II for details), to get 



A (q, Id) 



4:7ieuXr](q) 1 / sinhgd 



1 + 27re 2 M(0) 2q \cosh qd - 1 



(40) 
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where A (q, Id) and r)(q) are the two dimensional Fourier transforms of A (f,ld) and rj(f) 
respectively. In Eq. (40), M(0) is given by 

V ; J (2tt) 2 2g \coshqd- 1 / V ; 

and has the physical significance of an irreducible polarizability . Knowing the electric 
potential A (q, Id) (which is linear in u) we can integrate out the electrons and the harmonic 
electric potential fluctuations in the action (Eq. (22)) to second order in u. The result for 
the extra action per layer is 



Ske = J dt(u\y 



1 - 



2ne 2 M(0) 
1 + 27re 2 M(0) 



(42) 



where the first term in the square bracket is the large unscreened contribution calculated 
earlier (Eq. (31)), and the cancelling second term gives the reduction because of screening. 
Combining both these terms, we find m* to be 

m * = ft + 2,1,(0)) <43) 

We see from Eq. (43) that there is a 'dielectric' screening of the vortex mass, i.e. the 
factor in the denominator is a finite large number. This is due to the discrete level spectrum 
in the core, in contrast to the continuum of states in a metal or the near continuum in 
conventional superconductors giving 'metallic' screening. The core dielectric constant (e core ) 
can be evaluated using the variational wave functions of Eqs. (28) and (30). The screening 
reduces mj by a factor of about 50. The large dielectric constant can be understood as 
being due to the high polarizability of the core quasiparticle system. Approximately, the 
dielectric constant (e core ) is dimensionally given by the ratio of the core Coulomb energy and 
the excitation energy (e + _ = 2ei/2) necessary to create the particle-hole excitation which 
contributes to the polarization process leading to screening (see Appendix II). Thus, 

ecore ^ (^Coulomb/^excitation) - {(e 2 /0A+-} - 75 (44). 

The detailed calculation yields a dielectric constant of 53. Thus the screened or effective 
inertial mass of a vortex is rather small, being equal to 

m* = (mo/53) ~ 0.5m e (45) 

This is our main result, for the mass per layer, at T = 0, when the field is along the c-axis. 

We note that the mass Eq. (45) can be roughly estimated qualitatively by using the 
physical estimates Eqs. (32) and (44), i.e. 

m* ~ (h 2 k 2 F /e 1/2 ) /{(e 2 /e) /2e 1/2 } = {2k 2 F a G t)m e (46). 

This gives a value m* ~ 1.3m e close to the the result Eq. (45) of the detailed calculation! It is 
clear from the formal expressions, Eq. (29) and Eq. (43) for the effective mass, as well as the 
approximate form Eq. (46) that m* is effectively, the ratio of two polarizabilities, current- 
current and charge- charge. Vortex motion (in the rest frame of the vortex) gives rise to both 
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supercurrent fluctuations and electric potential fluctuations. They are connected because 
of gauge invariance. Since both polarizabilities involve the same energy denominator, this 
drops out of m* (Eq. (46)). It is also clear, from the occurrence of hp, a>o and £ in this 
equation, that the vortex mass is related to the carrier density, core size as well as Coulomb 
interactions. The small vortex mass in the cuprates is a consequence of the small core size 
and the low carrier density in comparision with conventional superconductors. 

In the next, concluding section, we discuss the approximations involved in the result 
obtained by us, its generalization, and the question of when vortex mass related effects may 
be observed. 

IV. Discussion and Conclusion 

A. Discussion of the Vortex Mass The microscopic calculation above uses an approx- 
imate variational order parameter A and corresponding Bogoliubov-deGennes amplitudes 
(u^Vfj) (Eq. (30)). An obvious improvement would be to solve exactly for these quantities 
given only the parameters A ,£ and the BCS relation £ = (Hvp/ttAq). This has been done 
numerically (22). There is of course a fair amount of uncertainty in these parameters for 
cuprate superconductors, quite apart from the question of whether an s-wave, BCS like or- 
der parameter with a conventional kinetic energy functional is at all appropriate for cuprate 
superconductors. However, within the s-wave BCS model, the exact expression Eq. (29) 
and (43) for m* can be evaluated, once u^, and their first spatial derivative are known 
for pi — 1/2. We have ignored the contribution due to transitions from the bound states 
to the continuum, satisfying the selection rule A/i = ±1. The reasons are the smallness 
of the matrix elements, the largeness of the energy denominator and very good screening. 
A rough estimate shows that these change the mass estimated by about 10-20%. We have 
also ignored contributions to the mass from polarisation processes involving the collective 
excitations of the superconducting state. In the case of a neutral Fermi superfluid it has 
been found by Niu et. al. (40) that the inclusion of these excitations, which correspond to 
long wavelength density fluctuations, leads to a finite vortex mass in contrast to the loga- 
rithmically divergent result obtained in our approach (see Section I.). However for the case 
of a charged superconductor, which is the primary concern of this paper, the corresponding 
mode is pushed up to the plasma frequency which is much larger than the other energies in 
the problem and is therefore unlikely to contribute to the vortex mass in a significant way. 

With increasing temperature, the gap A(T) and the inverse coherence length £(T) _1 
decrease. The structure of the vortex core also changes. In principle, one can repeat the 
T = calculation with temperature dependent input parameters. At low temperatures 
{ksT <C €i/ 2 ) this would roughly have the effect of increasing the effective inertial mass of 
the vortex as a function of temperature (Eq. 46). However, at higher temperatures (when 
€1/2 is of order or less than ksT) an additional contribution to the vortex mass will accrue 
from transitions involving thermally excited quasiparticles in the vortex core. A qualitatively 
new effect which arises in this regime is that due to quasiparticles scattering off the moving 
vortex, there is a damping of vortex motion. This dissipative term is generally included 
phenomenologically (it is linear in vortex velocity u and contributes an imaginary term 
to the action), though microscopic theories have been developed (25, 37). It is not clear 
whether the two effects viz. thermal renormalization of the effective mass, and dissipation, 
both due to thermal quasiparticles, are completely independent. Also, the regime where 
the bound core level spacing ei/ 2 is of order or less than /c B T is clearly very different from 
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the low temperature regime ei/ 2 >> ksT. We have not considered the former dissipation 
dominated 'high temperature' limit. 

There is the related question of adiabaticity, which is contained in the assumption that 
the order parameter of the system with one moving vortex is the same as that of the static 
vortex, but with r — > (r — ut). It is expected that the vortex motion would distort the gap 
function from its form in the static case. The distortions induced have to be determined self- 
consistently at every order of u. Simanek (26) has recently considered this question and has 
pointed out that the vortex velocity needs to be small enough such that (Uukp) < \ e i/2\ = 
(Ag/cp) or u < (uBcs)(A /ep) where ubcs is the BCS critical velocity i.e. TiUBcsk-F — A . 
The violation of this criterion means that the vortex motion can significantly distort the gap 
function as well as the magnetic field and supercurrent distribution. The changes induced 
by the vortex motion have to be self-consistently determined and it's effects included in the 
calculation of the vortex mass. However, at low vortex velocities (u <C ubcs), these effects 
are small and have therefore been ignored. 

We have made a Galilean transformation to the rest frame of the vortex to simplify the 
calculation of change in action to second order in u. Inclusion of the effects of the periodic 
lattice and impurity scattering is also possible within our formalism. Their interaction po- 
tential with the electrons becomes explicitly time- dependent and leads to inelastic scattering 
(38) of electrons whose effects have to be included while evaluating the correlation functions 
of Eqs. (24), (35) and (36). 

So far our discussion has been restricted to a single vortex. However, at larger magnetic 
fields, in the presence of a vortex lattice, additional contributions coming from the vortex 
lattice will have to be included. The periodicity of the pair potential broadens the localised 
quasiparticle levels into energy bands and new contributions to the vortex mass as well as 
the forces experienced by vortices are expected to arise from collective effects arising from 
the vortex lattice (15). A totally different approach becomes necessary in the presence of 
strong magnetic fields in the dense vortex limit near H C 2. The strong amplitude fluctuations 
which allow the dissociation of a Cooper pair make the dynamics of the order parameter 
diffusive in this regime. A calculation of the vortex mass will proceed from the Abrikosov 
solution (36) of the GL equations for a triangular vortex lattice. Both the vortex effective 
mass as well as the nature of the forces that drive the vortex dynamics in this regime are 
subjects that require further study. 

The core contribution to the vortex mass was estimated earlier microscopically by by 
T. Hsu (42). He obtained an answer which is of the same order of magnitude as the core 
contribution calculated by us in the abscence of Coulomb screening for parameters appro- 
priate to the cuprates. Hsu used the Bogulibov De-Gennes formalism to obtain the vortex 
acceleration in response to a transport supercurrent. The vortex mass is then deduced 
from this equation by a comparision with a hypothetical force equation obtained by setting 
the unknown vortex mass times the acceleration of the vortex equal to a Magnus force of 
the size suggested by Nozieres and Vinen using arguements of fluid hydrodynamics. While 
many aspects of the formalism are similar to ours we believe this method is unreliable for 
determining the vortex mass. The Coulomb screening effects which are found to be vital for 
giving rise to a small mass have been ignored. Even for the case of a neutral superfluid, the 
uncertainty in the size as well as sign of the Magnus force cast doubts about a procedure 
which relies on a knowledge of the Magnus force to deduce the core contribution to the 
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vortex mass. 

Recently, after completion of this work(15), a paper by Simanek (27) appeared where 
a calculation of the vortex mass using discrete core states is presented. This calculation is 
based on a TDGL functional, in terms of a single complex superconducting order parameter, 
which is derived in the presence of a moving vortex and leads to a mass which is dimensionally 
the same as our estimate for the unscreened mass. However, the all important screening 
effect which reduces m* by a factor of 50-75, has not been considered at all in Simanek's 
calculation. Further extensions (41) along the lines of Ref. (27) have appeared during the 
reviewing process. However, once again the substantial reduction in the vortex mass because 
of Coulomb screening has been ignored. 

The assumption of an s-wave like order parameter is not realistic. There is increasing 
direct experimental evidence for a strongly anisotropic order parameter, or an order param- 
eter with vanishing amplitude at some points in k (or r space) (28,29). Consider for example 
a d x 2_ y 2 like order parameter. The pair amplitude is nonlocal i.e. < ^|(r')^|(r) > vanishes 
for r = r', and has a dependence on the direction of (r — r') with a square symmetry. It 
is thus clear that the Bogoliubov-deGennes equations are nonlocal mixing different angular 
momentum eigenstates. Recent work by Volovik (30) and Ren and co-workers (35) (see also 
Ref. 44 and references therein) has suggested that vortices in superconductors with d x 2_ y 2 
symmetry have a non-zero s-wave component in the core of the vortex which vanishes at the 
vortex centre. Thus, the gapless bound state spectrum, which might have been expected 
for d-wave superconductors with lines of nodes in the gap function (in k-space) is absent. 
However, gapless excitations are available in the far region, where the s-wave component 
vanishes and are likely to give rise to strong dissipative effects so that the nature of vortex 
dynamics would be qualitatively very different. This is an area which needs much further 
work (see for example ref. 30 and 35). 

B. Observability of effects due to inertial mass 

It is clear that if there is no dissipation and no Magnus force, both of which produce a 
term in the action linear in velocity (10,11), the small inertial mass of a vortex would give rise 
to strong quantum effects. The vortices are bosonic particles whose degeneracy temperature 
is (eH/m*c)(h/k B )- This is of order 20K for an external field of 10T, if the effective mass 
m* is about 0.5m e ^. If this limit is realized, then several novel possibilities would arise, 
especially in strongly layered cuprate superconductors such as 2212. In these systems the 
vortex liquid phase extends to very low temperatures, especially in high magnetic fields (31). 
This liquid instead of becoming a solid, could on cooling become a quantum Bose liquid and 
then a genuine vortex superfluid (32). Such a vortex superfluid is a new ground state with 
unusual properties, most likely a new kind of insulator. The vortex superfluid could persist 
till T = 0, or freeze into a quantum solid, whose spectrum of collective excitations (phonons) 
would depend on the mass m*. The dynamics of vortices in this regime would be that of 
interacting bosons in a random potential. 

It is not clear however that the quantum Bose regime of the many vortex system is 
experimentally realizable. Firstly, at least for higher temperatures, there is strong dissipation 
which dominates the dynamics in both the quantum and classical regimes. Secondly there 
is a large Magnus force (9,10,11). If only the former were present, the mass could still be 
relevant for phenomena like quantum creep. If only the Magnus force were present, as is 
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believed to be the case in the cuprates where the onset of a dissipationless regime has been 
reported (5), the system of vortices is like that of bosons in a strong 'magnetic' field, the 
Magnus force being the analogue of the Lorentz force. The Hamiltonian of the system of 2d 
vortices can be written as 

H vort = E (P< " «,)72m* + (1/2) £ V{n - r 3 ) 

where (a x ,a y ) = (nhn/2)(—Xi,yi) and Vifi — fj) is the interaction between vortices. The 
a term is due to the Magnus force with n being the electron density per layer and m* is 
the vortex effective mass. The 'magnetic' field associated with this Magnus term is rather 
large, the cyclotron frequency being about 0.7 eV for m* ~ 0.5m e . Thus the Landau level 
separation is large, and the vortex system is in the lowest Landau level with a low filling 
fraction of (n v /n) where n v is the vortex density and n the electron density. The magnetic 
length of the system, i.e. the cyclotron orbit size is rather small ~ 7 A so that the dynamics 
is that of the guiding centre; the inertial mass is irrelevant. Even in the strong Magnus force 
limit, a large vortex mass has been shown (10) to give rise to quantum effects in phenomena 
involving vortex tunneling. In particular the semi-classical action develops a linear in T 
dependence, which would reflect in the observed rate of flux creep at low temperatures. In 
the language of our paper, a large vortex mass would result in a reduction of the cyclotron 
frequency, mixing in higher Landau levels and thus enhancing quantum effects. However, 
the rather small value of the vortex mass obtained by us implies that this scenario is actually 
not realised. The main uncertainty in vortex dynamics is the actual size of the Magnus force. 
The contribution of the bound states i.e. of localised quasiparticles to the Magnus force and 
the effect of disorder on it are major unsettled issues; there are a number of suggestions 
(33,34) that these could reduce, cancel or reverse the Magnus force! The dynamics of an 
isolated vortex with inertial mass, in the presence of a large Magnus force and dissipation 
had been investigated recently (9-12). Another possibility is that additional Magnus like 
forces could arise from the pair potential in the dense vortex lattice limit (15) with an 
opposite sign. However, there is a lack of a clear microscopic theory. There is a growing 
body of experimental evidence based on quantum creep (3), Hall measurements (6) and a.c. 
electromagnetic response (14) that the Magnus force is actually much smaller than current 
theoretical estimates (9, 10, 11). In that case, there exists the intriguing possibilty of the 
formation of a correlated quantum Hall fluid of the bosonic vortices at low temperatures 
(39). With the mounting evidence in the cuprates for a superconducting order parameter 
which has d x 2_ y 2 symmetry and quasiparticles whose mean free paths could be very long for 
T < T c , a realistic picture of this whole field awaits a microscopic calculation of the inertial 
mass, the Magnus force on a moving vortex and dissipation of its momentum for a d-wave 
superconductor at low temperatures. 
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Appendix I 

In this appendix we outline the evaluation of the correlation function (Eq. 24) which de- 
termines the unscreened core contribution to the vortex mass. Using the Bogulibov transfor- 
mation (Eq. 25) the field operators in Eq. (24) can be rewritten in terms of the quasiparticle 
operators (7^) to give 

J dtj dfj dr'D(f,r';t) (All) 



ITLn 



where the correlation function D(f,r';t) is given by 

i dx A ' i dx 



where 

pim = *m 7 j(*)7,(07i(oh„(o)]> (Ai.2b) 

. Here fl(r) = u^(r) and f^(r) = v*(f) and the summation with respect to pi, v, A and rj 
runs over both positive and negative values. The correlation function in Eq. (AI.2) is easily 
evaluated using the diagonalized Hamiltonian (Eq. (26)) to yield the expression given in 
Eq. (27). 

We will now derive the selection rule mentioned in Section III. To find mj from the 
correlation function of Eq. (27) we need to integrate with respect to the co-ordinates f and 
r'. The r -integration requires the evaluation of the integral 

Substituting the explicit forms of u^(r) and u ll i(f) (Eq. (28)) we find that 

h = ^J drrf-(r) J d<f> cos 

+h(pi' - 1/2) J drf-(r)f-,(r) J rf0sin0e i( ^'^ (AIA) 

The angular integrals in Eq. (AI.4) are zero unless pi — pi' = ±1. This together with the 
constraint pi > 0, pi' < implies that the only non-zero contribution to rriQ comes from 
= —pi' = 1/2. Making use of this selection rule, we find, after a little algebra, that the 
expression for m* reduces to Eq. (29). In arriving at this relation, we have used Eqs. (26b) 
and (26c). The only remaining task is to evaluate the matrix element g x occuring in Eq. 
(29). Substituting the explicit functional forms (Eq. (30)) into Eq. (29b) we find that 

% S^dxxe^MkF^l-MkF^x)^ + k F jJ' (k F {jx)} 
9x 2£ J™dxxe-*iJZ(k F Zx) + Ji(krfx)] 

Evaluating the dimensionless integrals on the R.H.S. of Eq. (AI.5) we find that for param- 
eters appropriate to the cuprates (k F £ ~ 4.47 ) 

I g x |^ 1.12| (AI.Q) 

The largest contribution to g x in Eq. (AI.5) comes from the term involving kp^J'^kp^x). 
However, unlike conventional superconductors, the relative smallness of the dimensionless 
parameter implies that the other term cannot be ignored. 
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Appendix II 

In the following, we present details of the calculation of the large reduction in the core 
contribution to the vortex mass due to Coulomb screening. To solve the Poisson equation 
(Eqs. (33) and (34)), it is necessary to determine the source term on the R.H.S. of Eq. (34) 
by integrating the correlation function xt £ (^-, r ') over the co-ordinate r'. Making use of the 
Bogulibov transformation ( Eq. (25)) and the diagonalised Hamiltonian (Eq. (26)) we find 
(after some algebra) that 



unocc. occ. 



ox/- "a o ^Kv)Uv{f) ,^.hdv*(r') „,-> s hduJr')i , ,att^ 
X ^(r,r') = 2 E E 77 _ \ M^)--^ ~ K^j-^j-L] + h.c. (AII.l) 

The operator k-j^ behaves like cos0 (sin0) and as before, the integration over the co- 
ordinate r' (see Eq. (AI.4) and the discussion following it) gives the selection rule fi — u = ±1 
which together with the constraint yt, > and v < implies /i = —v = 1/2. We therefore 
find 

/ ox- i u i -2^ /2 (r)^ /2 (r) r ( ^h du 1/2 (?) 

/ Xi {r,r)dr'= - / dr'v 1/2 (r')- — h h.c. [All. 2) 

In writing Eq. (AII.2) we have used Eqs. (26b) and (26c). Now substituting the explicit 
variational forms for Ui/ 2 (r) and Vi/ 2 (r) (Eqs. (28) and (30)) we arrive at Eqs. (37) and 
(38). To proceed further, we have to solve the Poisson equation and get a self-consistent 
solution for the scalar potential A (r,z) of the form assumed in Eq. (39). We will now 
consider the other term on the R.H.S. of Eq. (34). This term represents the screening 
charge induced by the Coulomb potential consequent to the electron density change induced 
by the vortex motion; it has to be determined by integrating the product of x°°(r, r')and 
A (r',ld) with respect to the co-ordinate r'. Since the latter has an angular dependence of 
the form sin we once again find that the only process which contributes to the polarization 
involves a transition from the highest occupied state (v = —1/2) to the lowest unoccupied 
state (pL — 1/2). We thus find, 

J d?x?{?, r')A (r>, Id) = J dr'A (r', ld)r)(r>) (AII.3) 

Combining Eqs. (37) and (AII.3) with the Poisson equation (Eqs. (33) and (34)) we get 

VM ° (r '^ = erj(f) V 5(z - ld)[u\ + & - [ dr'rj(r')A (r' : Id)] (AIIA) 

4:71 l 2 J 

Transforming to Fourier space, this can be rewritten in terms of the corresponding Fourier 
components as 

[-q 2 - k 2 ]A (q, k) /^v^ ■? ;ji e 2 i|(g] ^ f dq' -> 2ixm 



4tt 1 2d (2tt) 

(All.h) 

This is an integral equation for the scalar potential Aq. To solve for Aq, we find it convenient 
to introduce the quantity 
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Substituting Eq. (AII.6) in (AII.5) and making use of the property X(k) — X(k— ^p) for 
any integer to, we solve for X(k) to get 

_ -AneuXM(k) J2i e~ ikld 

X{k) ~ l + 2ire*M(k) {AnJ) 

where 

M(k) = [ JtlMl ( ^ d ) {AIL8) 



(2tt) 2 2q \ cosh qd — cos kd / 
Substituting Eqs. (AII.7) and (AII.8) in Eq. (AII.5) we can now solve for A (q, k) to get 

Mi, k) = ? + k2l+ ™> M{k) £ exp { -,U d ) (.4/7.9) 

Fourier transforming the above equation with respect to the wave vector k we get A (q, Id) 
(Eq. (40)). 

We now consider the action Eq. (22). We first consider the electric field energy. Inte- 
grating by parts, and making use of the Poisson equation (Eq. (33)) we get 

JdzJ d ^ VA ^ z))2 = ^E/ dr'A (f,ld)[u X nr,r')-e X ^(f,r')A (r',ld)} 

(All. 10) 

On integrating out the electrons to second order in the vortex velocity u and the Coulomb 
potential A we get the effective action 

/f -* — ?y 2 -» p 2 -» -* -» 

dt J drdr>[— X r(?,r') - -A (f,ld) X T(r,r')A (r',ld) + ueA (f,ld) X ^ x (r,r')} 

(AII.ll) 

Combining Eqs. (AII.10) and (AII.ll) we find that the effective action for the system, 
to second order in u, is given by 

S KE = T,J dt J dfdP[^fxi x (r, r>) + yA (f, ld) X T{r, ?)\ (AII.ll) 

The first term in Ske is the unscreened core contribution to the vortex mass evaluated earlier 
(Eq. (24)) while the other term represents the reduction because of Coulomb screening. We 
now substitute the explicit forms for xf x an d X°i x ■ Doing a calculation very similar to the 
one leading to Eqs. (37), (38) and (AII.2) we find 

J dfdPxF(r,P) = -2A 2 (All. 13) 



Substituting Eqs. (37) and (AIL 13) into Eq. (All. 12) we get 

euX 

t I dt[u*X* + 

i 



Ske = Y,J dt[u 2 \ 2 + ^ E / df v (r)A (f, Id)] (AII.U) 



Fourier transforming the second term in Eq. (AIL 14) and substituting the expression for 
A (Eq. (40)) we finally arrive at Eq. (42). 
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The only remaining task is to evaluate the polarizability M(0). Using Eq. (38a), we find 



that 



\ e V2 



d rf 1/2 (r)fy 2 (r) sin0e 



—iq-r 



(AII.15) 



On substituting the variational forms (Eq. 30) for fy 2 ( r ) an d /i/2( r ) an d integrating over 
the angular co-ordinate this reduces to 



2 27ri g y 



/■OO 

/ drrJ Q {k F r)J 1 {k F r)J 1 {qr)e- rli (All. 16) 

j0 



^ e V2 ^1/2 9 

Combining Eq. (AIL 16) with the expression for M(0) (Eq. (41)) we finally get 

An f , sinh(qd) 



M(0) 



A t/2 e i/2 J cosh(grf) - 1 



1%) 



where 



r (q) = / drrJ Q {k F r)J 1 {k F r)J 1 {qr)e- r/ ^ 
Jo 



(AIL 17) 



(A//. 18) 



Using Eqs. (30d), (43), (AII.17) and (All. 18) we are now in a position to calculate the 
'core dielectric constant' e core . We find that e core = 1 + 27re 2 M(0) is given by 



1 + 



where 



sinh( xd/t;) , 2 , . 

-/>) 



and 



roo 

L i = 4 / dx - , . , , , . 

jo cosh(xct/4j — 1 

/•OO 

L 2 = / dxxe~ x Ul(k F ^x) + J±(k F £x)) 
Jo 

roo 

f(x) = / dyyJ (k F £y)Ji(k F £y)J 1 (xy)e~ y 
Jo 



(All. 19) 

(AII.20) 
(All. 21) 

(AII.22) 



On evaluating these expressions (Eqs. (AII.19), (AII.20) and (AII.21) numerically we find 



^core — ^3. 
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Figure Captions 

Fig. 1 Feynman diagrams that contribute to the phase-only action functional of Eq. 
1) at T = in the clean limit on integrating out the electrons. The thick lines correspond 
to the fermions, the wavy lines indicate the density fluctuations p = (hO — 2eA ) and the 
dashed line stands for the supercurrent fluctuations j = (hS/9 — 2eA/c). 

Fig. 2 Additional diagrams that contribute to the Galilean invariant action functional 
of Eq. 7) on integrating out the electrons. The thick lines correspond to the fermions, the 
wavy lines indicate the density fluctuations p = (h9 — 2eAo) and the dashed line stands for 
the supercurrent fluctuations j = {KS/6 — 2eA/c). 

Fig. 3 Spectrum of bound states within the core (within the range — A to A ) and of 
continuum states outside it, with the angular momentum quantum numbers. The levels are 
appropriate to the parameters mentioned in the text for a cuprate superconductor at T = 0. 
The transition allowed by selection rules is shown by an arrow. 
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